Bound states in the continuum driven by AC fields 
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We report the formation of bound states in the continuum driven by AC fields. This system 
consists of a quantum ring connected to two leads. An AC side-gate voltage controls the interference 
pattern of the electrons passing through the system. We model the system by two sites in parallel 
connected to two semi-infinite lattices. The energy of these sites change harmonically with time. 
We obtain the transmission probability and the local density of states at the ring sites as a function 
of the parameters that define the system. The transmission probability displays a Fano profile 
when the energy of the incoming electron matches the driving frequency. Correspondingly, the local 
density of states presents a narrow peak that approaches a S function in the weak coupling limit. 
We attribute these features to the presence of bound states in the continuum. 
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I. INTRODUCTION 

At the dawn of quantum mechanics, von Neumann and 
Wigner constructed a spatially oscillating attractive po- 
tential that supported a bound state above the potential 
barrier X This truly localized (square integrable) solution 
of the time-independent Schrddinger equation is referred 
to as a bound state in the continuum (BIC). Much later, 
Stillingcr and Herrick reexamined and extended von Neu- 
mann and Wigner ideas ^ They analyzed a double excited 
atom model, where BICs were formed and live forever de- 
spite the interaction between electrons. They arrived at 
the conclusion that BICs may be a physically realizable 
phenomenon in real atomic and molecular systems. In 
this context, Friedrich and Wintgen discussed in a sys- 
tem of coupled Coulombic channels and, in particular, a 
hydrogen atom in a uniform magnetic fields These au- 
thors interpreted the formation of BICs as the result of 
the interference between resonances of different channels. 
The discovery of a BIC induced by the interaction be- 
tween two particles in close proximity to an impurity has 
been recently reported by Zhang et al., 4 where the state 
can be tuned in and out of the continuum continuously. 

The advent of nanofabrication techniques has made it 
possible to devise and fabricate quantum devices whose 
electronic properties are similar to those of atoms and 
molecules. When the size of the device is comparable to 
the de Broglie wavelength, one or more degrees of free- 
dom are quantized and the electron wave function is spa- 
tially confined. The similarity to atomic systems paved 
the way to experimentally validate the existence of BICs 
in artificial nanostructures. Capasso et al. measured 
the absorption spectrum at low temperature of a GalnAs 
quantum well with Bragg reflector barriers produced by 
a AlInAs/GalnAs superlattice^ A well defined line at 
360 meV in the spectrum was attributed to electron ex- 
citations from the ground state of the quantum well to a 
localized level well above the AlInAs band edge. Never- 
theless, this state cannot be regarded as a true BIC but 
a bound state above the barrier since it is a defect mode 



residing in the minigap of the superlattice, as pointed 
out by Plotnik et al& Thus, the experimental validation 
of the existence of truly BICs in quantum systems still 
remains a challenging task. Furthermore, it is worth to 
mention in this context that the analogy between pho- 
tonic systems in the paraxial regime and quantum sys- 
tems has facilitated the study£~— and subsequent experi- 
mental observation of BICsj* 

Electronic transport in mesoscopic and nanoscopic sys- 
tems can be also influenced by the occurrence of BICs. 
Ndckel investigated theoretically the ballistic transport 
across a quantum dot in a weak magnetic fields Reso- 
nances in the transmission were found to grow narrower 
with decreasing the magnetic field, and eventually they 
become BICs as the magnetic field vanishes. Fabry-Perot 
interference of quasibound states of two open quantum 
dots connected by a long wire cause the occurrence of 
BlCsiii The resulting state is nonlocal, in the sense that 
the electron is trapped in both quantum dots at the same 
time. Interestingly, controlling the size of one of the 
quantum dots makes the electron flows or get trapped 
inside the dots. BICs in parallel double quantum dot 
systems were found to be robust even if electron-electron 
interaction is taken into account^ Recently, Gonzalez et 
al. have demonstrated that not only quantum dots based 
on semiconductor materials but also on graphene can 
support BICsJ^ All these features stimulate the interest 
of BICs to develop new applications in nanoelectronics. 

In this work we extend the notion of BIC to the domain 
of time-dependent potentials. Our aim is twofold. First, 
we explore the possibility of the occurrence of BICs when 
the quantum system is driven by a time-harmonic poten- 
tial. Second, we introduce a physical realizable system 
which opens a novel possibility to reveal the existence 
of these exotic states in transport experiments. As a 
major result, we show that the transmission and corre- 
spondingly the conductance at low temperature signals 
the occurrence of BICs as Fano resonances. When the 
system is driven by an AC field the BICs survive and the 
existence is revealed by dynamic Fano resonances in the 
transmission probability. Remarkably, it turns out that 
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their energy can be tuned by changing the frequency of 
the field. Therefore, the conductance at low temperature 
presents a minimum when the BIC crosses the Fermi level 
by varying the driving frequency. 



II. QUANTUM RING UNDER AN AC 
SIDE-GATE VOLTAGE 



The system under consideration is a two dimensional 
gas of noninteracting electrons in a quantum ring, shown 
schematically in Fig. QJi). The ring is connected to two 
leads (source and drain). A side-gate voltage V±(t) 
breaks the symmetry of the upper and lower arms of 
the ring and acts as an additional parameter for con- 
trolling the electric current, as recently suggested for 
graphene-based nanorings . 14 ' 15 We assume that the side- 
gate voltage can be modulated harmonically in time with 
frequency w. 
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FIG. 1. a) Schematic diagram of the quantum ring with a 
side-gate voltage V±(t) connected to source (S) and drain (D). 
b) Equivalent lattice model with two time-dependent site en- 
ergies e± (t) at sites labeled 0± and two other sites with index 
±1 attached to semi-infinite chains. 



In order to study electron transport across the quan- 
tum ring, we mapped it onto a much simpler yet non- 
trivial lattice model, depicted in Fig. [Ud). We replace 
the actual quantum ring by four sites of a lattice within 
the tight-binding approximation. Two sites (0±) have 
time-dependent energies e±(t) and the other two sites, 
labeled ±1, are connected to semi-infinite chains. Time- 
dependent site energies are given by s+(t) = ±2AcoscJt. 
To avoid the profusion of free parameters, we assume a 
uniform transfer integral and vanishing site energies ex- 
cept at sites 0± , without loosing generality. The common 
value of the transfer integral will be set as the unit of en- 
ergy and we take h = 1 throughout the paper. 



A. Time-independent side-gate voltage 

To gain insight into the possible occurrence of BICs 
in the system, we consider the time-independent case by 
setting lj — for the moment. An incoming plane wave 
ij)™{t) = exp[i(kj— Et)], with energy E = —2 cos A: within 
the bands of the leads, will be partially transmitted in 
the form ipf(t) = t exp{i(kj — Et)). The lattice period 
is set as the length unit. It is a matter of simple algebra 
to obtain the transmission amplitude in this case 
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The transmission probability T(E) = |io| 2 presents a dip 
around the band center and vanishes at E = 0. 

In the weak coupling limit, namely A — > 0, the trans- 
mission probability shows a Fano profile 1 ^ close to the 
band center, T(E) = E 2 / (E 2 + A 4 ). The width of the 
dip scales as ~ A 2 . Poles of the transmission ampli- 
tude have a simple physical interpretation as the natural 
eigenstates of the scattering potential. — The poles occur 
at a complex energy whose real part gives the energy of 
the state and the imaginary part is related with its de- 
cay rate. There are four poles in the case under study. 
Two poles correspond to defect modes in the gap and, 
consequently, cannot be identified with BICs. However, 
the other two poles reside at the band center with an 
imaginary part equal to ztiA 2 when A< 1. In the weak 
coupling limit A — > the poles correspond to truly bound 
states. Therefore, the Fano resonance of the transmission 
amplitude (fTJ) signals the occurrence of BICs at the band 
center. 

To get a better understanding of the nature of the BICs 
at the band center we also calculate the local density of 
states (LDOS) at sites 0±, pa(E). Close to the band 
center, the LDOS is proportional to \ip + \ 2 + \ip~ | 2 , where 
is the wave function amplitude at those sites. The 
LDOS is given approximately as 
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The first term is nothing but the transmission probabil- 
ity, vanishing at the band center. However, the second 
term approaches Ait8{E) in the limit A — > 0, indicating 
the existence of a truly bound state with energy E = 
located at sites 0±. 



B. Time-dependent side-gate voltage 

We now turn to our main goal, the occurrence of 
BICs when the side-gate voltage depends harmonically 
on time. The time-dependent Schrodinger equation for 
the amplitudes ipj(t) reads 



(3) 
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where the index runs over the nearest- neighbor sites 
of j and the dot indicates the derivative with respect to 
time. Using the Floquet formalism, the solution can be 
expressed in the form 
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where E n = E + no; and n is the sideband channel index. 
Since we are interested in electron transmission across 
the ring, we take the following ansatz for the coefficients 
A n j in the expansion (|4"aj) 

{ 5 n0 e ik ^ +r n e~ ik ^ , j < —1 , 
A n ,i = { ft, 3 = 0± , (4b) 

[ t n e iknj , j > 1 • 

Inserting this ansatz in ([3]) leads to the dispersion re- 
lation E n — —2 cos fc„ where k n is real if E n lies within 
the band, i.e. \E + nuj\ < 2. In addition, we obtain 
t n = ft + fn = r n + ^nOi ensuring current conservation. 
Finally, one also gets 

A 2 A 2 

tn-2 ~ — t n +2 = 4id n0 sinfeo , (5a) 



where for brevity we define 
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The continued fraction approach developed in Ref. [l8| 
allows us to obtain numerically the contribution of all 
channels to the transmission. But if the coupling of the 
ring to the AC side-gate voltage is weak (A < 1), only 
the lowest order sidebands are significant. Then we keep 
five channels and assume that t n vanishes if \n\ > 3. 
Equation ([5]) implies that t±\ = in this approximation 
and 
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The transmission amplitude in the elastic channel is 
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to = 4i sin fco 
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Once the transmission amplitudes have been calcu- 
lated, we can obtain the transmission probability from 
the general expression 
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where the sum runs over the propagating channels, 
namely those channels for which E n — E + nuj lies within 
the band of the leads. Assuming that the sidebands are 
active for conduction, the transmission probability reads 

A 4 sin ko 
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A J 



\a 2 \ 2 Ef sinfeo 
sin fc_2 



|a_2| 2 £' 2 1 sinfeo 



(8) 



III. RESULTS 

The five channels approximation discussed above pro- 
vides a closed analytical expression for the transmission 
probability (|5J|. Admittedly the resulting expression is 
still involved and must be evaluated numerically for the 
chosen parameters. We will present some simpler expres- 
sions latter, valid in the weak coupling limit, but for the 
moment we are interested in the closed expression ©. 
We assume that coupling A is smaller or in the order 
of the frequency u> hereafter. Figure [5] shows the results 
when the frequency is cj = 0.5 and A = 0.1 (recall that 
the amplitude is actually twice this value). The pro- 
nounced dip observed at the band center in the static 
case (u = 0), commented above and shown in the figure, 
is absent if w ^ 0. Similar results can be obtained for 
other set of A and u parameters. In fact, transmission 
at E = is unity and the quantum ring becomes trans- 
parent at this energy. But the most salient feature of 
the transmission when the side-gate voltage oscillates is 
the occurrence of two symmetric and narrow dips, at en- 
ergies close to ±w. Remarkably, the transmission never 
vanishes in the range of energy plotted in Fig. [5] and at 
the dips only drops at about 0.5. Actually transmission 
probability vanishes but only at the band edges, as occurs 
in the static case too. The inset shows an enlarged view 
of one of the dips for two different values of the coupling 
A. It is quite apparent that the minimum transmission 
is slightly smaller than 0.5 and it is reached at an energy 
close but not exactly equal tow. In fact, transmission at 
E = u is exactly equal to 0.5 for any value of A. 




FIG. 2. Transmission probability as a function of energy at 
oj — (red dashed line) and ui = 0.5 (blue solid line) for 
A = 0.1. The inset shows an enlarged view of the transmission 
probability when E ~ cu at A = 0.01 (green dashed line) and 
A = 0.1 (blue solid line) for oj = 0.5. 



It is important to mention that we also solved numer- 
ically the general equation ([5]) to obtain the transmis- 
sion probability from ([7]), increasing the number of the 
sideband channels. The results were the same to those 
obtained within the five channels approximation, when 
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A is not too large. Thus, we can confidently use this 
approximation in our analysis. 

The transmission probability is an even function of en- 
ergy. Therefore, for concreteness we now focus in the 
energy region close +uj, when A is small. If uj is not 
large, we can take sinfc±2 — sinfco ~ 1 in (jHJ. In ad- 
dition, the term containing a.2 in Eqs. f|6b[) and © is 
negligible under these assumptions. After lengthly but 
straightforward algebra, the transmission probability re- 
duces to 
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Since uj and E — uj are not large, this expression can be 
further approximated by the following Fano profile 
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To facilitate the comparison with the static case, we have 
introduced an effective coupling A c ff in such a way that 
A 2 S = £±(i) = A 2 /2, where the bar indicates the average 
over one time period. Notice that the Fano factor q = 
—uj/A becomes independent of the coupling A. Figure [3] 
compares the exact results obtained from ©, for two 
different values of the coupling A and to = 0.5, and the 
Fano profile (|9b[) . In spite of the simplicity of the Fano 
profile and the assumptions we made, the agreement is 
remarkable in both cases. 
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FIG. 3. Transmission probability as a function of the param- 
eter x defined in (I9bjl . for two values of A and uj = 0.5. The 
Fano profile is a good approximation to both curves. 

In contrast to the static case, transmission at the dips 
remains finite when the side-gate voltage is harmonically 
modulated in time (see Fig. [5]). Therefore, it is not clear 
at this stage whether the dips are due to the occurrence of 
BICs in the system. To answer this question we consider 
again the LDOS at sites 0±. After time-averaging over 
one time period, one gets 
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Using the five channels approximation, the LDOS close 
to right dip (E ~ uj) becomes 
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and a similar expression is obtained for the left dip (E ~ 
—uj), replacing uj by — uj in (|10j) . Therefore, in the weak 
coupling limit A — > 0, the LDOS reduces to poHJ(E) ~ 
S(E + uj) + S(E — uj). In analogy with the static case, 
we claim that the two singular peaks in the LDOS at 
energies ±w are due to a new type of BICs arising by the 
interaction with the AC field. 



IV. CONCLUSIONS 

In summary, we have introduced and studied a novel 
type of BICs in systems whose energy levels are mod- 
ulated harmonically in time. To be specific, we have 
considered a quantum ring subjected to a side-gate volt- 
age oscillating in time with frequency uj and studied its 
transport properties in a fully coherent regime. We come 
to the important conclusion that the BICs supported by 
the quantum ring in the static case survive under har- 
monic modulation of the side-gate voltage. The two BICs 
driven by the AC field have energies ±uj and they re- 
veal themselves in the transmission, and consequently in 
the low-temperature conductance, as dynamic Fano res- 
onances. The position of the BICs inside the spectral 
band can be continuously tuned by varying the driving 
frequency and eventually they could be expelled out of 
the continuum when uj is larger than 2 in units of the 
transfer integral, i.e. when they approaches the band 
edge. Similar control of the energy of static BICs have 
been recently demonstrated by adding weak nonlinearity 
to semi-infinite systems 1 ^ or by varying the interaction 
between particles^ However, besides the different origin 
of the BICs, our proposal seems to be more advantageous 
for the experimental validation of these exotic states. In 
this regard, it is still an open question to what extend 
electron-electron interactions would mask the effect in 
a real experiment. Ztiko et al. have shown that the so- 
called dark states in parallel double quantum dot systems 
are robust against interactions within a Hubbard model, 
at least in the Kondo regime*^ These states correspond 
to the BICs of our present work in the static case, which 
makes us confident to expect also the BICs driven by AC 
fields to be robust in an interacting system. 
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